Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 30 August 2009 (MN KT^ style file v2.2) 



Lorentz Factor Constraint from the very early external shock of the 
gamma-ray burst ejecta 

; Yuan-Chuan Zou^'^, and Tsvi Piran^ * 

^The Racah Institute of Physics, Hebrew University, Jerusalem 91904, Israel 

^School of Physics, Huazhong University of Science and Technology, Wuhan 430074, China 

30 August 2009 



ABSTRACT 

While it is generally agreed that the emitting regions in Gamma-Ray Bursts (GRBs) move 
ultra relativistically towards the observer, different estimates of the initial Lorentz factors, 
To, lead to different, at times conflicting estimates. We show here that the quiet periods in 
which the signals goes down below the instrumental thresholds, put strong upper limits on 
the values of Fq. According to the standard internal-external shocks model an external shock 
should develop during the prompt stage. This external shock radiates in the hard X-rays to 
soft gamma-rays bands and this emission should be seen as a smooth background signal. The 
observed deep minima indicate that this contribution is negligible. This limits, in turn, Fq. We 
obtain upper limits on Fq for several bursts with typical values around hundreds. We compare 
these values with those obtained by the other methods, which typically yield lower limits. The 
results are marginally consistent leaving only a narrow range of allowed values for Fq. 
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1 INTRODUCTION 

The combination of fast variability and a non-thermal spectrum 
lead to the well known "compactness problem" and indicates that 
the emitti ng region of GR Bs move relativistically towards the 
observers (iRudermanI 1 1 975h (or that they are extremely nearby). 
With the realization, in the early nineties, that GRBs are cos- 
mological and the development of the fireball model, relativistic 
motion became an essential ingredient of GRB modeling. Rel- 
ativistic motion was confirmed, first by observations of scintil- 
lations at the radio afterglow of GRB 970508 (Goodman 19^; 
IWaxman. Kulkarni & Fr ail 1998) , later by afterglow modeling 
( Sari, ^iran & NaravanI 11998.: .vanParadijs Kouveliotou & WijersI 



bOOd : Penaitescu & Kumarl l2002h and finally in direct observa- 
tions of superluminal motion in the afterglow of GRB 030329 
(iTavlor et al. 2004; Oren, Nakar & Piran 2004). 

These observations were all done during the afterglow phase 
in which the ejected material has already been slowed down by its 
interaction with the surrounding material and the typical Lorentz 
factor was of order 5 or less. However, the compactness problem 
indicates that the initial Lorentz factor, Fq, with which the jet was 
ejected from the inner engine and with which it was moving during 
the prompt GRB phase, was much larger. 

In spite of this progress Fq, which is crucial to understand 
the underline physics both of the inner engine and of the emission 
process is still unknown. Several methods have been proposed to 
estimate it. The methods vary in complexity and in robustness and 
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depend on different details of the "internal-external" shocks model 
and/or on the development of an forward-reverse shocks system 
during the early phase of the afterglow: 

• The simplest, most direct and robust method uses the com- 
pactness problem: The optical depth for high energy photons (^ 
GeV) to escape from the emitting region without being annihilated 
by softer 7-rays (sub-MeV) should be les s than unity. This leads to 
a lower limit on the initial Lorentz factor dFenimore. Epstein & Hoi 
I993I: |PiranllT995l: IWoods & Loebl 1 19951 : tearing & hardindll99l 
Lithwick & Sarill200lh . It can be used to set an exact value for the 
Lorentz factor if an upper energy cuto ff is observed. However, this 
was not seen so far. Using this method jLith wick & S^(l200lh ob- 
tained lower limits of Fq (typically of order of several hundred) for 
several bursts. 

• Within the external shocks scenario the peak of the 
light curve corresponds to the decelerating time of the eject a 
(iMeszaros & Reeslll997l : ISari & Piraniri999 J bl: iKobavashi ll200Qh . 
Observations of this peak provide, therefore, a way to estimate Fq. 
Rykoff et al. ( 2009) constrained Fo using this method for several 
bursts who se prompt optical emission was observed by ROTSE- 
III. Zh ang et al.l (12006) argued that the deceleration time should be 
prior to the shallow decay phase, and used this to obtain lower lim- 
its ^100 for several bursts. 

• Comparison of the early X-ray and optical emission (that 
arise from a reverse-forward shock system) enables us to estimate 
the Lorentz factor hy fitti ng parameters of the emitting regions 
(ISari & Piranlll999al : lzhan g. Kobavashi & Meszaros 2003). As an 
example, using this method iSari & PiranI (11999al) constrained Fo of 
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GRB 990123 to be - 200. iMolinaril (l2007h and lJin&FanI (l2007h 
estimated using the early optical afterglow, To, to be ^ 400 for 
both G RB 060418 and GRB 060607A. Recently. Eue. Fan & Weil 
(l2009h used the reverse-forward shock model to determine the ini- 
tial Lorentz factor for several well observed bursts and obtained 
typical values of 300. 

• Thermal emission escapes from the fireball's photospheric 
when it becomes optically thin. The observations of a thermal com- 
ponent w ould allows us, therefore, to infer the Lore ntz factor at 
this time (iNakar. Piran & Sarill2005l : IPe'er et al.ll2007h . Assuming 
that such a compo nent was indeed observed and using this method 
iPe'er et al estimated To to be hundreds . 

In this work we propose yet another method to estimate Fq. 
We work within the "internal-external" shocks scenario accord- 
ing to which internal shocks that arise due to collisions of the 
ejected shells within the relativistic outflow produce the observed 
prompt gamma-rays, while the interaction of the merged shells 
with the surrounding medium produces the afterglow (see IPiranI 
I200I for a review). However, an external shock begins to develop 
even during the prompt phase. It is caused by the outermost shell 
that sweeps up the external medium. This external shock radiates 
and produces an underlying smooth component. Usually, this extra 
component will be sub-dominant when the internal shocks emis- 
sion is strong. However, as this emission is not observed even 
during troughs of the light curves it must be weaker than the 
detection threshold in cases when a deep minimum in the light 
curve is observed. For example. Swift BAT ha s a sensitivity limit 
of ^ 10~ ergscm" s~ ( Gehrels et al ]|20^ (corresponding to 
/zy,iim ~ 4 X 10~^^ergscm~^ Hz~^). We show here that the 
strength the early forward shock emission depends sensitively on 
the Lorentz factor at this stage (which is roughly the initial one). 
Using this we set upper limits on the initial Lorentz factor for vari- 
ant bursts. 

We derive the radiation flux density from the forward shock, 
and the resulting constraint on Fo in section |2] We examine other 
methods for constraining Fo in section [3] We calculate the con- 
straints for several selected individual bursts using three methods 
in section]?] and discuss the implications of the results in section |5] 



2 MODEL 

Our model is based on a basic ingredient of the "internal-external 
shocks model". Within this model, while internal shocks are going 
on producing the prompt gamma-rays the outermost shell that is at 
the front of the ejecta begins interacting with the surr ounding mat- 
ter and an a reverse-forward shocks system develops (Sari & PiranI 
I1995I : IZhang. Kobavashi & Meszarosl l2003l : iNakar & PiranI l2004l) . 
The reverse shock propagates back into the front of the ejecta and 
the external shock propagates into the surrounding matter. We con- 
sider here the emission from this very early reverse-forward shocks 
system. The contribution of the forward shock should appear in soft 
gamma-rays as a smooth and continuous emission with an increas- 
ing signal. However, in many cases the observed signal decreases to 
very low values, which can be even below the detection limit of the 
observing instrument. We use this to constrain the initial Lorentz 
factor. Depending on the environment, we consider two cases: in- 
terstellar medium (ISM) and a wind. 



2.1 An ISM 

The reverse shock can be relativistic (RRS) or Newtonian (NRS) 
depending on a density ratio between the ejecta and th e surrounding 
matter. The relevant case depends on the parameter JSari & PiranI 
.1995. : Nakar & Piran 2004): 
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where / = {?>E / Ai^nrripC^Y^^ is the Sedov length, Fo is the initial 
Lorentz factor of the merged shell, which is also the Lorentz factor 
of the prompt gamma-ray emitting region, Ao is the initial width 
of the shell, whic h is related to t he duration of the pulse by ~ 
Ao/(2c) 0.1s dWu et al.l2003h . E is the isotropic kinetic energy, 
n is the matter number density, and rup is the proton rest mass. The 
conventional notation Q = x 10^ is used throughout this paper. 
^ < 1 leads to RRS, while ^ > 1 corresponds to NRS. 

Typical values of the parameters have been used in eq. ([TJ. 
Unless ?7 or A are very large, <J ^ 1, so we only consider the 
NRS case. As the Newtonian reverse shock radiates near the op- 
tical band, the main contribution to the X-ray to soft gamma-ray 
emission arises from the forward shock, which we examine now. 

The deceleration radius for the outermost shell is 



Rd = 
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Consequently, the deceleration time is t 
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^s, which is much longer than the time 



of the first minimum in which we are interested. So we can con- 
sider the scaling-laws for the very early external shock before it has 
been decelerated. During this phase the shell coasts with a constant 
Lorentz factor Fo, collecting the medium to radiate X-rays and 
the shock accelerated electrons are cooled mainly via synchrotron 

emission. 

Following ISari. Piran & N aravan ( 1 998h (see also in e.g. 



iNakar & Piraif l2004l : IZou^ Wu & Dai.,2005ir we calculate the dy- 
namics and radiation from the forward shock. The internal energy 
density of the forward shocked material is e = ^rfmpC^n 
60no Fo,2 ergcm"^, the magnetic field in the comoving frame 



h B — y^STTCBe ^ 12nQ ^0,2 Gauss, where cb is the 

equipartition factor for the magnetic energy density. The peak 
spectral power is ,max — 

4.7 X 10"^^ (1 + z)nj e| _iFg,2ergHz-^s"\ where ar is 
the Thomson cross section, qe is the electron charge. The peak 



observed flux density is then 
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where A^e is the total number of emitting electrons, and D is 
the luminosity distance. The synchrotron cooling Lorentz fac- 



tor IS 7c 
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z)t'^^ (1 + y)~'^ £3^-1 (where tco is the comoving time 
scale and Y is the Compton parameter for synchrotron self- 
Compton scattering). This corresponds to the cooling fre- 



(1 + 



-1 Tpj^qeB 
27rmeC 



quency i^c 

-3 _2 

trons is 7m ^ _^ee_ 
typical synchrotron frequency is 
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3|3Y, and p is the index of power law distributed electrons. The 
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Synchrotron- self absorption frequency is Va ^ 2.8 x 10^ (1 + 

13 9 8 6 28 

z)~~ n§ t| ^B,-i ^0% fo^" i^a < i^c < i^m and Ua ~ 4.6 X 
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10^ (1 + z) 
which are always well below Vm and Vc. 

Collecting the above expressions, the observed flux density fv 
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for p = 2.43 and > (z^a, z^c, z^m)- One can see that fi, de- 
pends very sensitively on the initial Lorentz factor To . This is par- 
tially because of the number of emitting electrons A^e, for A^e = 
47r/3i?^n oc V%t%n, behaves like V%. 

For different parameters, especially for different values of To, 
V may be either greater or smaller than Vm or Vc. Figure [T] shows 
rough spectra and the relations between the different frequencies 
for 3 different values of Fq. Even though the observed frequency is 
not always larger than Vc and Vm, it is never significantly smaller. 
Therefore, equation ([3]) is roughly acceptable for the whole range 
of parameters of interest. For simplicity, we use only this equation 
in the following discussion. 

If the observed flux density immediately after the first GRB 
pulse is /i/,c, The emission of the early external shock should not 
exceed this value, requiring ^ fv,c- In other cases, there is 
no signal immediately after the first pulse, when < fu,iim, 
where fu,iim is the limiting flux density of the observing instru- 
ment. We use fu,iim to represent both quantities in the following. 




Figure 1. The sketched synchrotron spectra of the forward shock emission 
before the deceleration. The lines indicate for Fq = 400, 200 and 100 from 
top to down respectively. 



2.2 A wind environment 

For a wind environnient, th e distinction between RRS and NRS 
depends on (IWu et al.ll2003h 



Using 



10 erg cm ^ Hz ^ , this inequality gives a 



strong constraint on Fo : 
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where the Sedov length in a wind case is / = Eo/{47vAmpC^), 

Fo < 340(1 + ^) V^^Hm,-28^2l<*e^! 1 es*i^^2f ^0^(1 + ^) ^ (^F ^^^^ parameter satisfying n(r) = Ar'^, and A = 3 x 

' 10^^ cm~^. For typical parameters the NRS case is also most 



Note that the limit practically depends rather weakly on all the other 
parameters. For completeness, we also present the expression for a 
genrael p: 



likely. The deceleration radius is 
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and the corresponding deceleration time is t^^d = 30(1 + 
2;)^53r^2A~^^sec, which is also longer than the duration of a 
typical first pulse. The observed flux density of the forward shock 
before deceleration is: 



[3(p-2)/(p-i)]-(^-^)(i + y)]- 
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As the value of p ranges between 2 and 3 generally, and as the 
overall expression depends on 2ff4^ ^o is weakly dependent on 
the other parameters for a general p. 
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^ Note that for simplicity, the above relations are given for p = 
2.5. For other values of p, the indices will change accordingly 
but not too much. For example, if p = 2.2, the scaling law 



in the case u > {iym,i^c,i^a) fsee lZou. Wu & Pail (l2005h for other 
cases). As the wind density is much higher than the ISM density, 
the flux density is much higher than the one expected for an ISM 
(given by eq. (|3])). For the typical parameter A^ = 0.1, the number 
density at a given radius R = 2Fgct© - 6 x lO^^Fg 2^0 cm IS 
about 10^ larger than the ISM, and consequently, the flux density 
is about 10^ stronger, as it is proportional to n^^^. For p : 
initial Lorentz factor can be expressed as 



2.5, the 
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in the gamma-ray band is fu,ma^{i^c/i^m)~^^~-^''^'^{i^/i^c)~^^'^ or 
fu,ina^ii^c/i^m)~^^~'^'^^'^ii^/i^c)~^^'^, which is lowcr with higher value 
of p. Therefore, the choice of p = 2.5 is conservative when evaluating the 
mini mal flux density, as the usual value of p is around 2.2 ( Achterberg et al. 
I2OO I). Furthermore, the value of p in the range of [2,3] doesn't change the 
emission much. See also eq. {5} for a general p. 



The limit is even stronger than the one obtained in ISM case. This 
is essentially because of the denser medium for the typical wind 
parameter. However, the density profile is uncertain at small values 
of r (that are relevant at this stage) and it is not clear if this limit 
is valid. In any case, this is a strong argument against very dense 
environment near GRBs. 
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3 OTHER CONSTRAINTS 



3.3 A photospheric thermal component 



As mentioned earlier several other methods to constrain the initial 
Lorentz factor have been suggested. We compare the constraints 
obtained using the different methods to show the consistency of the 
overall model. As most other methods give lower limits, combined 
with ours, we obtain a overall stronger constraint. 



3.1 Compactness and the optical depth for pair production 

Observations of high energy photons from a GRB imply that the 
emitting region is optically thin for pair production of these high 
energy photons with the lower energy sub-MeV gamma-rays. This 
leads to a limit (Lithwick & Sar i 2001) (denoted limit A) 
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where (3 is the photon spectral index {fE ^cm ^erg ^), 
^max is the maximal energy of the observed high energy photons, 
and T is defined as 

. _ (ll/180)crT/:>'(mec')-^+V 
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(11) 



where D is the luminosity distance, is the variability time scale 
of the prompt emission, and /i = /(MeV)~^^^s~^cm~^ is the 
number of photons per second per square centimeter per MeV at 
energy of 1 MeV. 

When high energy observations are not available, we can still 
obtain a bound by assuming that the observed higher energy spec- 
trum at around an MeV can extend up to the energy where photons 
can just annihilate with themselves in the same spectral shape. This 
leads to another constraint dLithwick & Sari. 200 L) (denoted limit 
B): 
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Typically, limit A yields a tighter constraint than limit B. 



3.2 The deceleration time 

The deceleration time of the external forward shock, t^^2 corre- 
sponds to the peak in the afterglow light curve. It depends weakly 
on the overall energy and the external density and most sensitively 
on the initial Lorentz factor. At the deceleration time , the Lorentz 
factor i s a half of the initial Lorentz factor. This yields dSari & PiranI 
Il999bl) : 
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where is the deceleration time, that corresponds to the peak in 
the afterglow light curve. 

By employing an efficiency parameter fj = ^^,iso/^k,iso, 
where ^k,iso and ^7,iso are isotropic equivalent kinetic energy and 
the prompt emitted photon energy respectively, we get 



To ^ 220^;/L,3,r;:o\4'<'/'C/'(1 + z) 



x3/8 



(14) 



For a wind environ ment, the initial Lorentz factor is 
dPanaitescu & Kumai]l200Ql) 



where = 3 x lO^^cm"^ is the wind parameter. 



(15) 



A thermal component, should come from the fireball's photosphere 
(where the optical depth is unity). Detection of such a compo- 
nent would enable us to determine the initial Lorentz factor using 
(.Nakar. Pira n & Sari 2005; Pe'er et al. 2007): 
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where IZ ■ 
emission. 



{Fbb/cjT^u 



2mpC^n_ 
and Fbb is the thermal component 



3.4 The reverse-forward shocks system 

During the early afterglow, the predicted reverse-forward shock 
may give a direct clue to determine the initial Lorentz factor, as 
the Lorentz factor of the unshocked ejecta (that is involved with 
the reverse shock) is roughly equal to the initial Lorentz factor. By 
fitting the early afterglow to the predictions of the reverse-forward 
shock model, one obtains different parameters within the shocked 
region including the Lorentz factor. This is the final Lorentz factor 
after the merging during internal shocks, which may be regarded as 
an average value of the initially ejected sub- shells. 



We have listed several methods. However, not all are prac- 
tical. Specifically, it is difficult to dig out a thermal component 
within the total emis sion as nee ded for the thermal component 
method (see however [loos'). Similarly, the identification of 
a reverse-forward shocks system is limited and available only for a 
few bursts. In the following we will use only the methods in sec- 
tions OCLSOl A or LSOIB) and[32](SP99). 



4 CASE STUDY AND A COMPARISON OF THE 
DIFFERENT METHODS 

We examine some well observed bursts to see whether the differ- 
ent constrains on the initial Lorentz factor are consistent with each 
other. The new constraint described here is based on the observa- 
tion of the first pulse and the following minimum. For subsequent 
pulses the forward shock is complicated. Therefore, the constraint 
is only suitable for those bursts which exhibit a short first pulse 
followed by a deep trough. In the following we consider only such 
bursts. To use the three methods, we need the bursts have the prop- 
erties: (i) A GeV signals (for LSOIA). (ii) A clear first pulse and 
a subsequent minimum (for this work, ZP09), with not long du- 
ration, otherwise, the external shock may have decelerated during 
the first pulse, (iii) A peak in the optical afterglow light curve (for 
SP99). However, it is very rare satisfying all the three requirements. 
Then we reduce it to any burst who has two of the above preper- 
ties as the sample to examine the Fq. We use typical parameters: 
/i/,iim 4 X 10~^^ergscm~^ Hz~^ as the detection limit if 
it is a Swift burst; fj 0.3, n = lcm~^, p = 2.5, Ce — 0.3, and 
= 0.1 if they are not determined. 

GRB 990123 

GRB 990123 was located at z = 1.60 (i:>L = 3.7 x 10^^ cm) 
with an isotropic eq uivalent e nergy in 7-rays ^7,iso ^ 3.4 x 
10^^ ergs fe ulkarni etal .11 19991) . It has been extensively discussed. 
lLithwick& Saril (i200l]) obtained a lower limit of Fo ~ 150 using 
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LSOIA and 180 using LSOIB (assuming there is a single power law 
in the spectrum of the high energy band). 

The afterglow hght curve can be fitted both by a uni- 
form medium and a wind-type medium, with parameters 
n = 0.004cm"^, = 0.075, = 4 x 10"^, for uniform 
medium and = 0.06, Ce = 0.08, = 7 x 10"^ for 
wind-type medium respectively, while the underly ing physical 
parameters are somewhat uncertain (|Panaitescull2005). The prompt 
BATS E light curve shows a valley at ~ 12sec ( Galama et al.l 



BATSE sensitivity in 50-300 keV, is 0.2cm"2s"^ 



Il999l) . 

(http://glast.gsfc.nasa.gov/science/instruments/ta blel-2.htm l), 
corresponding to a flux density of fu^ 

lim ^ 1.25 X 

10~^^ergcm~^s~^Hz~^. Using those parameters and equa- 
tions © and 12]), we obtain an upper limit of To < 1200 for a 
uniform medium and 410 for win d-type medium by o ur method. 

From the optical light curve ( Galama et al ]| 19991) . the deceler- 
ating time should be less then 0.1 day (~ 8640 s). (Though there 
were also optical observations before , they are generally consid- 
ered as reverse-forward shock signals dSari & Piranlll999ah .) Using 
equations ilAi and ([15]), and assuming the efficiency ry = 0.3, we 
obtain lower limit of 100 for uniform medium and 130 for wind- 
type medium respectively by method SP99. 

These constraints are con sistent with t hose o btained by light 
curve fitting, e.g. To = 270 bv lKobavashil (|200Q|) and To = 1200 
bv lWang. Dai & Lul (l2000h . 



GRB 021004 

ORB 021004 was observed b y HETE-II, and located at redshift z — 
2.32 {Dl ^ 5.8 X 10^^ cm) jpox et al.l2003l) . Lacking high energy 
data, we cannot method LSOIA. To use method LSOIB, we need 
the spectrum of the prompt emission. However, the spectru m fitting 
was a single power-law with photon index (3 ^ 1.64 (La mb et al.l 
l2002h . An extended single power-law with /3 < 2 indicates the 
most energy is hidden in the higher energy band and there must 
be a spectral break and an unknown high energy spectral index. 
Therefore, we cannot carry out the estimation by method LSOIB 
neither. 

The duration of the first pulse is ^ lOsec 
Chttp://space.mit.edu/HETE/Bursts/GRB021004/). HETE's Detec- 
tion thresholds for the French Gamma Telescope (FREGATE) is ^ 
3 X 10"^ergcm~^s"^ (http://space.mit.edu/HETE/fregate.html), 
which corresponding to fu^Um 1.2 x 10~^^erg cm~^ Hz~^ 
at 100 keV. The afterglow emission was fitted well by wind- type 
environment, with parameters — 0.6,^^,52 — 10, Ce — es — 
0.1 (Li & Chevalier 2003). Using equation ([9]), we obtain the upper 
limit of To < 210. 

From the optical light curve, which wa s first observed 537sec 
after the trigger of the burst (IFox et al.ll2003h . the decelerating time 
of the external shock should be shorter than this time. Using equa- 
tion ([15]), we obtain (by method SP99) To > 80. 



GRB 040924 

GRB 040924 was a short burst with duration Tgo = 2.39=b0.24 s.It 
was l ocated at z = 0.858 {Dl ^ 1.7 x 10^^ cm) ('Wierse ma et all 
l2005h. The peak energy of the time-integrated spectrum obtained 
by Konus-wind was Ep — 67 =b 6 keV. Howeve r, the spec- 
tral indices were not available jGolenetskii et al.|[2004l) . Therefore, 
method LSOl can not be used. 

The 7-ray isotropic equivalent energy was ^ 1.5 x 10^^ 



erg dFan et al. l2005l). The first opt ical observations was taken at 
- lOOOsec JWiersema etaDl2005h . which yields an upper limit 
of dec elerating time. As the parameters are not firmly determined 
Jpan et al. 2005), we choose the typical parameters (with HETE- 
II limit: /^,Hm ^ 1.2 x 10"^'^erg cm"^ Hz"S"^). We find 
To < 490 by method ZP09, and To > 120 using SP99. 

GRB 050401 

GRB 050401 was located at z = 2.9 {D l ^ 7.6 x 10^^ cm). Th e 
duration of the first pulse was ^ 6sec (De Pasquale et al.l l2005l) . 
From Konus-Wind observation (Golenetskii et al. 2005a), the 2nd 
peak of the prompt pulses had Ep = 119 =b 26 keV, a = 0.83, 
(3 = 2.37 and peak flux 2.45 ± 0.12 x lO-^ergcm-^s"^ (in 
the 20 keV - 2 MeV energy range). Without a direct high energy 
(^GeV) detection, we use method LSOIB to get a lower limit. Us- 
ing /i - 0.27cm"^s"^MeV"^ and f - 1.7 x 10^°, we get 
To > 110. Using the BAT limit and other ty pical parameters, and 
ZP09 we obtain To < 590. The lower limit iRvkoff etaP (l2009h 
using method SP99 is To > 900. 



GRB 050801 

GRB 050801 was located at z = 1.56 {Dl ^ 3.6 x 10^^ cm) 
There was no observed emission at time 8sec dDe Pasquale et al.l 
l2007h . The peak flux was 1.7 ± O.lph cm"^s"^ in 15-350 keV 
at Isec, and the time averaged spectral index was 13 = 2.0 =b 0.2 
dSakamoto et al. l 2005) (we use this as the index at the peak time). 
Taking the duration of the first pulse to be 2sec, we get /i ~ 
0.03cm-2s-^MeV"\ f - 1.5 x 10^ and To > 80 using 
LSOIB. Giving the typical parameter values of n, ee, and es, we 
get To < 420 using ZP09. The lower limit using SP99 is To > 500 
(lRvkoffetal.ll2009h . 



GRB 050922c 



GRB 050922c had T90 - 4.5s (Sakamo to et"aDl2008h . It was lo- 
cated at z = 2.198 (Jakobsson etal 2005i) {Dl ^ 5.5 x 10^^ cm) 

with fluence / = 3.1 x 10"^erg cm"^ in 30-400 keV dCrew et al.l 
l2005h . which corresponds to Ei^o ^ 3.7 x 10^^ erg. As the s pec- 
tral po wer-law index was (3 — 1.55 ±0.07 (< 2) (Golenetskii et al.l 
l2005bh . we cannot use method LSOl. Taking T90 as the duration of 
the first pulse and the typical parameters, we get To < 550 using 
ZP09. A lower limit using method SP99 is To > 350 dRvkoff et al.l 
l2009h . 



GRB 060607a 

GRB 060607a was locat ed at z = 3.082 {Dl ^ 8.2 x 10^^ cm) 
dNvsewander et al ]|200i. With no high energy observation and a 
BAT photon spectral index (3 — 1.45 =b 0.07 we cannot use ni ethod 
LSOl. The first peak ended at 15 sec dZiaeepour et al.l2008b . The 
isotropic equivalent energy was ^iso ~ 1.1 x 10^^ erg. A clea r 
optical peak was observed at ~ 180sec dNvs ewander et al.ll2009h . 
With typical parameters, the upper limit using ZP09 is To < 490. 
The inferred initial Lorentz factor using SP99 is ^ 410. 



GRB 060614 

GRB 060614 was located at redshift z = 0.12^ {Dl = 1.8 x 
10^^ cm) with an isotropic equivalent energy E^^^^o 2.5 x 10^^ 
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ergs ^Mund ell et alJ bOOvl) . Ther e was an exponential cu t-off at 
^ 300 keV for the intense pulse (iGolenetskii et alj|2006ah . which 
means there are practically no photons at higher energy. So we can 
not u se method LSOl. T he duration of the first pulse was about 
Isec jGehrels et alJl2006h . From afterglow modelling, the physical 
parameters are Ek ^ 6 x 10^° erg s, te ^ 0.12, cb ^ 0.0002, 
and n = 0.04cm"^ dXu et al.ll2609l) . Taking the instrument limit 
of Swift BAT, we find (using ZP09) To < 530. The peak time of the 
afterglow was ~ 3 x 10^ sec (Mundell et al . 2007). Taking this as 
the decelerating time, and ry = 0.3, we find (using SP99) To ^ 35. 
This pretty low value may imply that the peak of the optical after- 
glow for GRB 060614 did not occur at the de celeration time but it 
arised due to energy injection as suggested bv lXu et aP (l2009l) . 



GRB 061007 

This burst was located at z = 1.26 ( Dl ^ 2.7 x 10^ ^ cm). The 
duration of the first pulse was ^ 3 sec dOhno et al.ll2009h . The peak 
time was at - 39.5sec, with Ep = 498148 keV, a = 0.53lo S9' 
= 2.6lt° o^ and the peak flux 1.95to\lleTgcm-^s-^ 
(IGolenetskii et al.ll2006bh . Taking the duration of the peak pulse 
to be 0.3sec(see http://gcn.gsfc.nasa.gov/notices_s/232683/BA/ for 
the prompt light curves), we get /i ~ 3.5 cm~^s~^MeV~^, 
f - 6.7 X 10^\ and To > 160 using LSOIB. 

T he micro-physics pa rameters being uncertain dMundell et al.l 

I2OOI [ Schadv et al.l l2008h . we choose the normal value Ce — 
0.3, es = 0.1 and n — lcm~^, which don't affect the re- 
sult much. Using the Swift BAT limit and equation lU, we find 
To ^ ^^^n~^^^e~^_^lj^e~B/J^ with ZP09. The optical afterglow 
peaked at ^ 39sec, and the isotropic equivale nt gamma-ray en- 
ergy was - 1.4 X 10^^ erg (iRvkoff et al.ll2009h . We obtain, using 
equation ([T4l), Fp ^ 640 fylQ^^n^^^^, which is consistent with 
iRvkoffet all (12009 ). 

The later two constraints are inconsistent, and there is no much 
space to tune the parameters and reach consistency. However, one 
have to recall that method SP99 is somewhat crude and it is not 
clear that the inconsistency of less than a factor of 2 is significant. 
Moreover, a possible explanation for this contradiction could be 
that the Lorentz factor of the first ejected shell is less than 460, 
while the followed other shells move faster than the first one. After 
they merged, produce a single shell moving with a higher Lorentz 
factor. This may be an evidence that the outermost shell is acceler- 
ated during the prompt phase by other shells. 



GRB 080319B 

GRB 08 03 19B, the naked ey e burst, was located at redshift z — 
0.937 ( Vrees wiik eTaDl2008l) . Its duration T90 was 57sec. The 
peak flux was Fp 2.26=b0.21 x 10"^erg cm~^s"^ and the peak 
of the vFi, spectrum was Ep ~ 675±22keV (i.e., Vp ~ 1.6 x 
10^°Hz, consequently - 1.4 x 10~^^erg cm"^Hz"^s"^). 
The photon indexes b elow and abo ve Ep are 0.8551q qJ4 and 
3.59^0 32 respectively jRacusin et al.l l2008). With luminosity dis- 
tance D - 1.9 X lO^^cm, GRB 080319B had a peak luminos- 
ity Lp ^ 9.67 X 10^^erg g~^ and an isotropic equivalent energy 
E^,iso ^ 1.32 X 10^^ erg (^Golenetskii etal.'2008). The peak pho- 
ton fluence is roughly ~ 20cm~^s~^ (Racusin et al. 2008). Taking 
the duration of pulse ^ 2sec, we get /i ^ 2.3 cm~^s~^MeV~^, 
f - 3.9 X 10^° and Fi,min 50 using method LSOIB. The dura- 
tion of the first pul se was about 3sec, an d the valley after the pulse 
was about 1 mJy jRacusin etal]l2008h . With = 0.01, = 



0.2, - 6 X 10"'^ (iRacusin et al.ll2008l), we find, using Z P09, 
Fo < 580. From the optical light curves (IPandev et al.ll2009h . one 
can estimate the deceleration time of the afterglow to be less than 
lOOsec. Setting the ef ficiency parameter f] — 0.3, and the wind 
parameter 0.01 dRacusin et al.ll2008h . we find, using SP99, 

Fo > 810. 



GRB 080916C 

This burst was located at z 4.35 (Dl ^ 1.2 x 10^^ cm). To 
carry out method LSOl, we focus on the 2nd pulse, which took 
place between 3.6 - 7. 7 sec. Two photons with energy > IGeV 
were detected during this pulse. For the soft gamma-ray band, 
the peak energy was ^ 1170 keV, the photon spectral index 
in the higher band was ^ 2 .21, and the peak flux density was 
3.5 X 10~^ cm"^s"^keV~^ jAbdo et al.ll2009h . and then the flux 
at 1 MeV was /i - 48 cm"2s"^MeV"\ We findf - 1.34x10^^ 
and Fq > 870 us ing LSOl A. This is consistent with the result of 
lAbdo et al.l (120091) . while Fo > 490 by method LSOIB. 

Consider the first p ulse of this burst, which was during the 
period 0.004 - 3.58sec (lAbdo et al.ll2009l) . The external forward 
shock should no t exceed the observed average gamma-ray flux 
- 6.9 cm-^s"^ (lAbdo et a l. 2009) in energy band 50 - 300 keV. 
As the photon index was ~ 0.58, the flux density at 300 keV was 
~ 1.6 X 10~^^ergscm~^ Hz~^ s~^, which can be taken as a con- 
servative fu,iim. With typical parameters and using eq. (|4l), we get, 
using ZP09, Fi < 1130. The optical afterglow showed that the de- 
celeration time was less than 6 x lO^s ( . Greiner et a l. 2009). With 
E^,iso ^ 8.8 X 10^^ ergs (lAbdo etal 1120091) . 



we get the initial 



Lorentz factor Fo > 90 using SP99. 



GRB 090328A 



This burst was located at z = 0.736 (ICenko et al.ll2009l) (Dl ^ 
1.4 X 10^^ cm). The spectrum from T0+3.1sec to T0+29.7sec was 
best fitted by a Band function with indices a = 0.93 ± 0.02 and 
/3 = 2.2 ± 0.1, and peak energy of ^peak = 653 ± 45 keV. The 
fluence in this time interval was 8.09 ± 0.10 x 10~^ergcm~^ in 
the 8-1000 keV band and 9.5 ± 1.0 x 10~^erg cm"^ in the 8keV- 
40MeV band. The isotropic equivalent energy in the 8keV-40MeV 
band was ^7,iso = 2.3 =b 0.2 x 10^^ ergs. The 1-sec peak photon 
flux measured starting from T0-i-23.5sec in th e 8-1000 keV band 
was 18.5±0.5phs-^ dRau et al.l200"9ah . which corresponds 

to /i - 2.7 phs"^ cm"^ MeV"\ and therefore f - 4.2 x 10^° 
(taking the variability time scale to be Isec). The Fermi Large Area 
Telescope (LAT) had detec ted this GRB with emission observed 
up to a few GeV jMcEnerv et al. 200 9). However, the ar rival time 
was very uncertain, even up to 900sec (ICutini et al.l2009l) . Here we 
assume the observed highest photon is 5 GeV, and it was in the 
same time interval of the prompt soft 7-rays. We get Fo ^ 320 
using LSOl A and Fo ^ 130 using LSOIB. The first pulse was 
about 4.2sec dRau et al.ll2009ah . By taking typical parameter, we 
find Fo < 540 using ZP09. As the early full optical light curve is 
not available, we cannot carry out the SP99 method for this burst. 



GRB 090424 

GRB 090424 was located at z = 0.544 dChornock etal.ll2009h 
(Dl = 9.6 X 10^^ cm). The peak energy of this burst was 
Ep = 177 ± 3 keV, with a = 0.9 ± 0.02 and /3 = 2.9 ± 0.1 
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The 0.128-sec peak photon flux mea- 
sured at 1.4sec in the 8-1000 keV band was 137 ± 5ph s~^cm~^ 
(pulse duration 0.3 s ) dConnaughton et alJ l2009h . We get 
/i - 1.6cm-2s-^MeV-\ f - 4.0 x 10^°, and To ^ 70 us- 
ing LSOIB. Without observations of high energy (GeV) photons , 
method LSOIA cannot be used. 

The fluence (8-1 000 keV) over the entire event was 5.2 ± 
0.1 X 10~^ergcm~^ dConnaughton et al.] |2009V corresponding to 
the E^^iso 4 X 10^^ erg. We take the vall ey at 6sec (the first 
major pulse was ^ 6sec with a few sub-pulses dConnaughton et al.] 
I2OO9I) ) to constraint the initial Lorentz factor using ZP09. Using 
typical parameters, we get To ^ 300. 

The optical temporal i ndex varied from 1.2 during ^ 
100— ~ lOOOsec IXin et al .12009 ). This indicate deceleration time 
should be less than ~ lOOsec. Using method SP99, we get the 
To ^ 310. 



GRB 090510 

GRB 090510 was classified as short burst with a duration 0.5sec 
jHove r sten et al. _ 2009). It was located at a redshift z 0.903 di 
0.003 jRauetal.ll2009bh . corresponding to a luminosity distance 
1.8 X 10^^ cm. The integrated spectrum was well fitted by a Band 
function with a = 0.80 ± 0.03, = 2.6 ± 0.3 and Ep = 4.4 ± 0.4 
MeV, and the 8 keV to 40 MeV fluence was ^ 3.0x10"^ ergcm"^ 
llGuiriec et al. 2009) , corresponding to the total isotropic equiva- 
lent energy 6.4 x 10^^ erg. The peak photon flux was 80 cm~^s~^, 
corresponding to /4.4 ^ 24 cm~^s~^MeV~^ (at 4.4 MeV) and 
/i '-^ 204 cm~^s~^MeV~^ (which is the extension from the /3 
slope), SOT ~ 5.2 X 10^^. There were > 10 photons with energy 
> 1 GeV during the prompt phase dOmodei et al.ll2009l) . Using this 
photons and taking the duration of the sub-pulses is 0.1 sec, the min- 
imal Lorentz factor satisfies To > 960 using LSOIA and To > 340 
with LSOIB. The minimal energy to annihilate the high energy 

photon is ^max,an = (LmeC^) V^max 26ri^;;^^ -^o^^^MeV, 

which is well above the break point of the soft-gamma spectrum. 
Therefore, the photon index 2.6 can be safely used. 

The first valley occurred at ~ 0.1 sec. Using equation ©, we 
find To < 620 using Z P09. UVOT found an optical peak at - 
600sec (iKuin et al.ll200^ . Taking this as the deceleration time of 
the external shock, and assuming n — lcm~^, 77 = 0.3, we obtain 
Lo - 180 using SP99. 



The limits on the initial Lorentz factor obtained using different 
methods are summarized in Table [T] and depicted graphically in 
fig.El For a few bursts (GRBs 050401, 050801, 061007, 080319B, 
090424, 090510) the limits are inconsistent. One may wonder if this 
inconsistency is problematic. First, we should realize the "initial 
Lorentz factor" for the different methods point to different objects. 
In the method SP99, it is the "final" Lorentz factor after all the sub- 
shell merged. In method LSOl it corresponds to the specific shell 
which produces the GeV photons. In our method, ZP09, it is the 
Lorentz factor of the first shell. It is possible that different objects 
have different Lorentz factors even for the same event. 

Moreover, the methods that depend on different assumptions, 
may not be that accurate. First, all three constraints assume the re- 
lation R — 2ro c^t to obtain the emission radius. Method SP99, 
assume no energy injection, and it depends on parameter such as 
the density n and the gamma ray efficiency fj. When using LSOl, 
Ep should be less than Emax,an, and the high energy spectrum 



Table 1. Comparison of the initial Lorentz factor constraint of different 
methods. The inconsistent cases are in boldface. 



Burst No. 


z 


SP99 


LSOl 


ZP09 


medium 


GRB 990123 


1.60 


^ 100 


^ 180* 


^ 1200 


uniform 






^ 130 




^410 


wind-type 


GRB 021004 


2.32 


^ 80 


a,b 


^ 210 


wind-type 


GRB 040924 


0.858 


^ 120 


a,c 


^490 


uniform 


GRB 050401 


2.9 


^ 900 


^ 110* 


^ 590 


uniform 


GRB 050801 


1.56 


^ 500 


^ 80* 


^420 


uniform 


GRB 050922C 


2.198 


^ 350 


a,b 


^ 550 


uniform 


GRB 060607A 


3.082 


~ 410 


a,b 


^490 


uniform 


GRB 060614 


0.125 


- 35 


a,d 


^ 530 


uniform 


GRB 061007 


1.26 


- 640 


^ 160* 


^480 


uniform 


GRB 080319B 


0.937 


^ 810 


^ 50* 


^ 580 


wind-type 


GRB 080916C 


4.35 


^ 90 


^ 880* 


^ 1130 


uniform 


GRB 090328A 


0.736 


e 


^ 320* 


^ 540 


uniform 


GRB 090424 


0.544 


^ 310 


^ 70* 


^ 300 


uniform 


GRB 090510 


0.903 


~ 180 


^ 960* 


^ 620 


uniform 



a. no high energy (~GeV) observations; b. photon index is less than 2; 
c. no photon indices; d. exponential cutoff; e. no optical data available. 
* limit A of LSOl; * limit B of LSOl. 



GRB 090510 




Figure 2. The allowed space (in grey) for Fq for different bursts by three 
different methods. For each burst, from bottom up the methods are SP99, 
LSOl and ZP09 respectively. For the limits A and limit B in method LSOl, 
only the higher value is used. 

should obey the Band function. Finally for ZP09, particle's power- 
law distribution and early formed external shock are assumed. 



5 CONCLUSION AND DISCUSSION 

We have considered the very early emission from the external 
shock, that begins to develope already during the prompt stage and 
emits soft 7-rays. Such emission was not identified so far. It is pos- 
sible that we have seen such early soft 7-ray emission from the 
external shock, either as a single smo oth light curve like the one ob- 
served in GRBs 911022 and 920216 (iFishman etalll 19941) . or as a 
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smooth pulse follow ing some additi onal pulses like GRB 50525 A 
(lBlustinetalJ[2006 ) and 0809 16A dZiaeepour et aDl2009l) . or that 
the signal was hidden among numerous pulses like GRBs 911106 
and 911127, 920221 (Fishman et al. 1994). However, one cannot 
confirm that these signals arose from external shocks. On the other 
hand, in many bursts a clear strong minima, reaching the sensitivity 
limit of the detector that follows the first pulse puts a very strong 
limits on this emission. These limits constrain the physical param- 
eters of these GRBs and in particular the initial Lorentz factor. In 
cases when a clear minima is seen after the first prompt pulse this 
leads to a strong upper limit, typically of order of a few hundreds 
on the Lorentz factor. The exact value depends rather weakly on the 
sensitivity of the observing instrument, the distance, the density of 
the environment, the duration of the first 7-ray pulse and on the 
micro-physical parameters. 

In view of the insensitivity of the constraint to various param- 
eters, it is rather robust, provided that it is applicable and that the 
external shocks model is relevant at this stage. For example the very 
early external shock results from the interaction with matter that is 
rather close to the progenitor star and the environment is rather un- 
certain. These considerations rule out a dense wind profile that 
extends all the way to small distance from the progenitor as this 
would produce a too strong early external shock signal. Addition- 
ally, it applies to the Lorentz factor of the outermost shell, which 
could be slow relative to subsequent shells that follow. 

A comparison with two other independent methods to estimate 
the Lorentz factor reveals that the three methods are inconsistent 
(by a factor of up to 1.5) in 5 out of 14 cases considered. This factor 
of 1.5 may not be significant taking into account the uncertainties in 
some of the methods. In method SP99, the Lorentz factor depends 
on the total kinetic energy and the environmental density, while the 
kinetic energy is very uncertain. In method LSOl, the Lorentz fac- 
tor depends sensitively on the spectral index (3 which is uncertain 
especially in the higher energy band. The assumption of a spectral 
single power-law in the higher energy, used in LSOIB, may also be 
invalid. Moreover, we should notice that the different methods ac- 
tually address different "initial Lorentz factors": the "final" Lorentz 
factor after all the sub-shells merged (in SP99), the Lorentz factor 
of the shell emitting the highest energy photons (in LSOl) and the 
Lorentz factor of the outmost shell (in ZP09). 

Additional bursts, and in particular additional bursts contain- 
ing GeV emission detected by Fermi for which the compactness 
problem is most efficiently utilized will enable us, hopefully in the 
near future, to confront the very early afterglow constraint with the 
lower limits obtained by the compactness problem. Consistency be- 
tween the two will confirm that we are on the right track towards 
a resolution of how GRBs work, while a significant contradiction 
will pose yet another puzzle. 



ACKNOWLEDGMENTS 

We thank R. Sari, and Y. Z. Fan for the helpful discussion. This 
work is supported by an ERC advanced research grant and by the 
center of excellence in High Energy Astrophysics funded by the 
Israel Science Foundation by the Schwartzmann chair (TP) and by 
the National Natural Science Foundation of China under the grant 
10703002 (fYCZ). 



REFERENCES 

Abdo A. A. et al., 2009, Science, 323, 1688 

Achterberg A., Gallant Y. A., Kirk J. G., Guthmann A. W., 2001, 

MNRAS, 328, 393 
Baring M. G., & Harding A. K., 1997, ApJ, 491, 663 
Berger E., et al., 2003, Nature, 426, 154 
Blustin A. J., et al., 2006, ApJ, 637, 901 
Campana S., et al., 2006, Nature, 442, 1008 
Cenko S. B., et al., 2009, GCN circular, 9053 
Chornock R., et al., 2009, GCN circular 9243 
Connaughton V., et al., GCN circular 9230 
Crew G., et al., 2005, GCN circular, 4021 
Cutini S., et al., 2009, GCN circular, 9077 
De Pasquale M., et al., 2005, MNRAS, 365, 1031 
De Pasquale M., et al., 2007, MNRAS, 377, 1638 
Fan Y. Z., Zhang B., Kaobayashi S., Meszaros R, 2005, ApJ, 628, 

867 

Fenimore E. E., Epstein R. I., & Ho C, 1993, A&AS, 97, 59 

Fishman G. J., et al., 1994, ApJ, 92, S229 

Fox D. W., et al., 2003, Nature, 422, 284 

Galama T. J., et al., 1999, Nature, 398, 394 

Gehrels N. et al., 2004, ApJ, 611, 1005 

Gehrels N. et al., 2006, Nature, 444, 1044 

Golenetskii S., et al., 2004, GCN, circualr, 2754 

Golenetskii S., et al., 2005a, GCN, circualr, 3179 

Golenetskii S., et al., 2005b, GCN, circualr, 4030 

Golenetskii S., et al., 2006a, GCN, circualr, 5264 

Golenetskii S., et al., 2006b, GCN, circualr, 5722 

Golenetskii S., et al., 2008, GCN, circualr, 7482 

Goodman J., 1986, ApJ, 308, L47 

Goodman J., 1997, New Astron., 2, 449 

Greiner J., et al., 2009, arXiv:0902.0761 

Guiriec S., et al., 2009, GCN circular 9336 

Hoversten E. A., et al., 2009, GCN circular 9331 

Jakobsson, P., Fynbo, J. P. U., Paraficz, D., Telting, J., Jensen, B. 

L., Hjorth, J., & Cern, J. M. C. 2005, GCN Circular 4029 
Jin Z. R, Fan Y. Z., 2007, MNRAS, 378, 1043 
Kobayashi S., 2000, ApJ, 545, 807 
Kuin N. P M., et al., 2009, GCN circular, 9342 
Kulkarni S. R., et al., 1999, Nature, 398, 389 
Penaitescu A., Kumar P, 2002, ApJ, 571, 779 
Lamb D., et al., 2001, GCN Circular 1600 
Li Z. Y, Chevalier R. A., 1999, ApJ, 526, 716 
Li Z. Y, Chevalier R. A., 2003, ApJ, 589, L69 
Lithwick Y, Sari R., 2001, ApJ, 555, 540 
Mangano V., et al., 2007, A&A, 470, 105 
McEnery J., et al., 2009, GCN circular, 9044 
Meszaros P, & Rees M., 1997, ApJ, 476, 232 
Molinari E., et al., 2007, A&A, 469, L13 
Mundell C. G., 2007, ApJ, 660, 489 
Nakar E., Piran T., 2004, MNRAS, 353, 647 
Nakar E., Piran T., Sari R., 2005, ApJ, 635, 516 
Nysewander M., Reichart D. E., Grain J. A., Foster A., Haislip J., 

Ivarsen K., Lacluyze A., & Trotter A., 2009, ApJ, 693, 1417 
Ohno M., et al., 2009, arXiv:08 12.3737 
Omodei N., et al., 2009, GCN circular, 9350 
Oren Y, Nakar E., Piran T., 2004, MNRAS, 353, L35 
Panaitescu A., 2005, ApJ, 363, 1409 
Panaitescu A., Kumar P, 2000, ApJ, 543, 66 
Pandey S. B., et al., 2009, arXiv: 0904. 1797 



Lorentz Factor Constraint 

Pe'er A., Ryde R, Wijers R. A. M. J., Meszaros R, Rees M. J., 

2007, ApJ, 664, LI 
Piran T., 1995, Proceedings of "Some Unsoved Problems 

in Astrophysics", Princeton, Eds. J. Bahcall and J. Osriker 

(arXiv:astro-ph/9507114) 
Piran T, 2005, Rev. Mod. Rhys., 76, 1143 
Racusin J. L., et al., 2008, Nature, 455, 183 
Rau A., et al., 2009a, GCN circular, 9057 

Rau A., McBreen S., Kruehler T., Greiner J., 2009, GCN circular, 
9353 

Ruderman, M., 1975, in 7th Texas Symposium on Relativistic As- 
trophysics, Dallas, Tx., 1974, Ann. N.Y. Acad. Sci. 262, 164 
Ryde R, 2005, ApJ, 625, L95 
Rykoff E. S. et al., 2009, arXiv: 0904. 0261^ 
Sakamoto T., et al., 2005, GCN Circular, 3730 
Sakamoto T., et al., 2008, ApJS, 175, 179 
Sari R., Piran T., 1995, ApJ, 455, L143 
Sari R., Piran T., 1999a, ApJ, 517, L109 
Sari R., Piran T., 1999b, ApJ, 520, 641 
Sari R., Piran T., Narayan R., 1998, ApJ, 497, L17 
Schady R, et al., 2008, AIPC, 1000, 200 

Taylor G. B., Prail D. A., Berger E., Kulkarni S. R., 2004, ApJ, 
609, LI 

van Paradijs J., Kouveliotou C, & Wijers R. A. M. J., 2000, 

ARA&A, 38, 379 
Vreeswijk R M., et al., 2008, GCN circular, 7444 
Wang X. Y., Dai Z. G., Lu T., 2000, MNRAS, 319, 1159 
Waxman E., Kulkarni S. R., Frail D. A., 1998, ApJ, 497, 288 
Wiersema K., et al., 2005, A&A, 481, 319 
Woods E., & Loeb A., 1995, ApJ, 453, 583 
Wu X. R, Dai Z. G., Huang Y. R, Lu T., 2003, MNRAS, 342, 1131 
Xin L. R, et al., 2009, GCN circular 9225 
Xu D., et al., 2009, ApJ, 696, 971 

Xue R. R., Pan Y. Z., Wei D. M., 2009, arXiv: 0902. 2613" 
Zhang B., Pan Y. Z., Dyks J., Kobayashi S., Meszaros P., Burrows 

D. N., Nousek J. A., Gehrels N., 2006, ApJ, 642, 354 
Zhang B., Kobayashi S., Meszaros R, 2003, ApJ, 595, 950 
Ziaeepour H., et al., 2008, MNRAS, 385, 453 
Ziaeepour H., et al., 2009, GCN circular, 8237 
Zou Y. C, Wu X. R, Dai Z. G., 2005, MNRAS, 363, 93 



